Recognition of stable distribution with Levy index alpha close to 2 
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We address the problem of recognizing alpha-stable Levy distribution with Levy index close 
to 2 from experimental data. We are interested in the case when the sample size of available 
data is not large, thus the power law asymptotics of the distribution is not clearly detectable, and 
the shape of empirical probability density function is close to a Gaussian. We propose a testing 
procedure combining a simple visual test based on empirical fourth moment with the Anderson- 
Darling and Jarque-Bera statistical tests and we check the efficiency of the method on simulated 
data. Furthermore, we apply our method to the analysis of turbulent plasma density and potential 
fluctuations measured in the stellarator type fusion device and demonstrate that the phenomenon 
of L-H transition occurring in this device is accompanied by the transition from Levy to Gaussian 
fluctuation statistics. 
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I. INTRODUCTION 

Stable distributions, also called alpha-stable, or Levy 
stable are ubiquitous in nature due to Generalized 
Central Limit Theorem. It says that the stable distri- 
butions, like the Gaussian one, attract distributions of 
sums of independent identically distributed random vari- 
ables ,2j. Due to this reason, Levy stable distributions 
naturally appear when evolution of a system or result of 
an experiment are determined by a sum of many random 
factors. An important parameter characterizing stable 
distribution is its index of stability, or the Levy index 
a, < a < 2. If a is strictly less than 2, the proba- 
bility density function (PDF) of the stable law exhibits 
slowly decaying power law asymptotic behavior of the 



form 



("heavy tail"), hence, the variance diverges. 



Due to this reason, the stable PDFs appear naturally in 
the description of random processes with large outliers, 
far from equilibrium. The limiting case a = 2 gives Gaus- 
sian distribution possessing fast decaying tails and finite 
variance. In general, each stable distribution is charac- 
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terized by four parameters: Levy index a, skewness pa- 
rameter /?, —1 < /3 < 1, scale parameter 7, and location 
parameter 5. For the Gaussian distribution j3 is irrele- 
vant, and it is characterized by standard deviation and 
mean. The comprehensive theory of alpha-stable Levy 
distributions is presented in monographs 

Some measurable quantities obey Levy statistics ex- 
actly, such as, e.g., gravitational field of masses dis- 
tributed randomly in space, hitting times for one- 
dimensional Brownian motion, and points of arrival in 
two-dimensional one @ ; interesting examples can be also 
found in [3|, [7[ • Levy statistics may also appear asymp- 
totically due to Generalized Central Limit Theorem, like, 
e.g., in non-Brownian random walks with jumps and/or 
waiting times obeying heavy-tailed distributions, see the 
reviews [cl-fiol]; we mention here such illuminating ex- 
amples as circulation of dollar bills [ll[ , Levy flights for 
light in fractal medium called Levy glass [l2[ , and Levy- 
like behavior of the marine vertebrates in response to 
patchy distribution of food resources Stably dis- 

tributed random noises are observed in such diverse ap- 
plications as plasma physics (density and electric field 
fluctuations [lj, [la]), stochastic climate dynamics 1161. 
physiology (heartbeats [r?|), electrical engineering |18| . 
biology [19( , and economics [2fJ . These and a lot of other 
observations of stably (or stable-like) distributed quan- 
tities require reliable methods of random data analysis 
which allow to detect stable distributions (or distribu- 
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tions which belong to their domains of attraction) and 
estimate their parameters. 

There are at least three procedures for estimating 
Levy stable law parameters: (i) the maximum likelihood 
method based on numerical approximation of the Levy 
stable likelihood function [2l|; (ii) the quantile method 
using tabulated quantiles of Levy stable laws (22[; and 
(iii) the method using regression on the sample charac- 
teristic function (23[. All presented methods work well 
assuming that the sample under consideration is indeed 
Levy stable. The regression method is both fast and ac- 
curate, thus we use it in our analysis below. However, 
if the data come from a different distribution, these pro- 
cedures may mislead more than the Hill |24f and direct 
log-log scale tail estimation method |25|, > which focus 
only on the tail behavior without assuming parametric 
form for the whole distribution function. 

The present paper, motivated by the analysis of plasma 
data (see Section V below), addresses the issue of recog- 
nizing stable distributions with the Levy index close to 
2. In this case, if the data sample is not large enough, 
the shape of empirical PDF is close to a Gaussian, and 
both log-log scale analysis and Hill estimator give over- 
estimated value of the Levy index for the number of ob- 
servation less than 10 6 [26|], see also examples in Figs. 
3-7 of Ref. (25[. In applications, especially in physics 
and biology, the number of observations is often less. In 
particular, we deal with plasma data containing a few 
thousand of data points. Thus, one has to be cautious 
in that domain of the Levy indices. Here we propose 
a certain testing algorithm which, to our belief, could 
be helpful for the data analysis in that case. It combines 
"visual inspection" of Levy-type behavior with a few sta- 
tistical tests, namely, the Anderson-Darling and Jarque- 
Bera. Being direct and relatively simple, they, however, 
are rarely used in physical and biological applications, 
to our knowledge. We demonstrate effectiveness of the 
proposed testing procedure by using simulated data with 
known Levy stable distribution. As an application, we 
analyze the PDFs of plasma fluctuations measured in the 
fusion plasma device and detect an interesting change of 
plasma fluctuation statistics from Levy stable to a Gaus- 
sian one in the so-called L-H transition which occurs in 
the plasma device during operation. 

The structure of the paper is as follows. In Sec. |TT] 
we first plot empirical PDFs for simulated data with the 
Levy index close to 2 and 5000 observations, and make 
sure that the difference from the Gaussian distribution is 
almost not visible. Then we plot the cumulative fourth 
moment of the simulated data as a function of obser- 
vation number and demonstrate that this simple visual 
test may help to distinguish between the Levy stable and 
Gaussian PDFs. In Sec. IHII we outline the algorithm used 
to distinguish between the Levy stable and Gaussian dis- 
tributions and demonstrate how this algorithm works for 
simulated data. This is accompanied by a short sketch of 
the statistical tests which we employ for the analysis. In 
Sec. HVI we apply the testing algorithm for the analysis of 



the experimental data from plasma physics. Finally, in 
Sec. |V]we add a few remarks on the two issues: the pro- 
posed testing procedure and change of statistics observed 
in Sec. [TV] 

II. SIMULATIONS 

In order to demonstrate behavior of the Levy stable 
distribution with Levy index close to 2 we first turn to 
the analysis of simulated samples. The procedures of sim- 
ulating stable random variables one can find in [26j j . 
In Fig. Q] we present in the log-log scale the right tail of 
the empirical PDF of the symmetric stable distribution 
with parameters a = 1.98, /? = 0, a = 1 and /i = 
together with the empirical PDF of the Gaussian distri- 
bution with mean equal to zero and variance equal to 
2. Let us mention, the Gaussian distribution with zero 
mean and variance equal to 2 corresponds to the Levy 
stable one with a = 2, a = 1 and fi = (the /3 param- 
eter is irrelevant in this case). In both examined sam- 
ples number of observations is set to 5000. Moreover, in 
Fig. Q]we also show the straight line with the slope pa- 
rameter —2.98 corresponding to asymptotics of the Levy 
stable PDF with a — 1.98. As we observe, the differ- 
ence between the two empirical PDFs is hardly visible. 
Therefore, we propose to use a simple visual test to dis- 
tinguish between the distributions. Namely, we calculate 
the empirical cumulative fourth moment (ECFM) of the 
simulated data sets, which for a sample of observations 
{xi, . . . , x n } is defined as follows: 

1 k 

C ^= fcE^-^) 4 ' k = l,2,...,n, (1) 

i=l 

where x is the mean of the random sample. 

The formula (JI]) can be calculated for any sample ob- 
tained from an arbitrary probability distribution. For 
a fixed k it forms a random variable. For distributions 
with finite fourth moment (e.g., Gaussian), the ECFM, 
as a function of k, converges to a constant, whereas for 
distributions with infinite fourth moment (e.g., Levy sta- 
ble with a < 2) it diverges to infinity. The latter, for a 
finite sample, can be observed as an irregular chaotic be- 
havior. Such reasoning stands for the idea of the first 
step of our algorithm. This simple test, which is pre- 
sented in Fig. [T] (inset), clearly indicates that there is 
a noticeable difference between Levy and Gaussian dis- 
tributions. For the former the ECFM does not tend to 
a constant value with number k of observations increas- 
ing, and behaves chaotically, while for the latter one for 
large numbers of observations it goes to the theoretical 
fourth moment that in this case is equal to 12. Now, if we 
take stably distributed variables with smaller Levy index, 
say, a — 1.9, and the same total number of observations 
n, then the difference between the empirical Levy stable 
PDF and Gaussian PDF becomes more visible; this is 
shown in Fig. [5J Of course, in this case we also observe 
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the chaotic behavior of the ECFM for stably distributed 
variables, see the inset in Fig. [2J In Fig. [3] we demon- 
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FIG. 1. (Color online) The two empirical PDFs calculated for 
simulated samples from Levy stable distribution with param- 
eters a = 1.98, ft = 0, a = 1 and fi — 0, and from Gaussian 
distribution with mean equal zero and variance equal 2, are 
shown in the log-log scale. Number of observations in both 
samples is 5000. The straight line shows the asymptotic slope 
of the Levy stable PDF with a — 1.98. Inset: Empirical 
cumulative fourth moment for the samples considered. 

strate the right tail of the empirical PDF (in the log-log 
scale) for totally skewed Levy distribution with a = 1.98, 
ft = 1, cr = 1, and fi = 0. Number of observations is the 
same as in Figs. [1] and [3J i.e., n — 5000. Similarly to 
the symmetric case /3 = 0, in Fig. [3] we also plot the 
PDF of the Gaussian distribution with mean equal zero 
and variance equal to 2, and the slope of the asymp- 
totics of the Levy stable PDF. The difference between 
the considered PDFs is almost not visible, similar to Fig. 
[TJ However, in contrast to Fig. [TJ the visual test based 
on ECFM does not indicate clearly a difference between 
the analyzed distributions, since for both empirical PDFs 
the ECFM tends to a constant, see the inset in Fig. |3J 
This example demonstrates that obviously there is a need 
to use a more advanced technique for recognition of the 
Levy stable distribution with a close to 2. 



III. PROCEDURE OF TESTING 

In this section we present an algorithm how to recog- 
nize the difference between the Gaussian distribution and 
the Levy stable distribution with the index a close to 2. 
It is based on the empirical fourth moment studied in de- 
tail in Section Ulland Anderson-Darling and Jarque-Bera 
statistical tests. 

Goodness-of-fit statistical tests, in general, are con- 
structed to check whether a hypothetical distribution can 
be rejected for the given data or not. Their idea basically 
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FIG. 2. (Color online) Empirical PDF calculated for simu- 
lated sample from Levy stable distribution with parameters 
a = 1.9, P = 0, a = 1 and /i = and for Gaussian dis- 
tribution (with mean equal zero and variance equal 2) in the 
log-log scale. Number of observations in both samples is 5000. 
The straight line shows the asymptotic slope of the Levy sta- 
ble PDF with a = 1.9. Inset: Empirical cumulative fourth 
moment for the samples considered. 



relies on checking how probable is reaching the value of 
some statistic calculated for the data if the underlying 
distribution was the hypothetical one. 

Statistics usually measure the distance between the 
empirical and the fitted analytical cumulative distribu- 
tion functions. For a sample of observations {x\, . . . , x n } 
the empirical cumulative distribution function F n (x) is 
a piecewise constant function starting from zero with 
jumps of size 1/n at points Xi. The distance is us ually 
measured either by a supremum or a quadratic norm [28| . 

The most well-know supremum statistic is the 
Kolmogorov-Smirnov (KS) statistic. It is just the supre- 
mum of the set of distances: 



D = sup\F n (x)-F{x) 



(2) 



where F(x) is the analytical cumulative distribution 
function. 

The Cramer- von Mises family of statistics incorporates 
the idea of the quadratic norm. It is defined by: 



n / {F n (x) - F(x)} 2 ip(x)dF(x), 



(3) 



where tp( x ) is a suitable function which gives weights to 
the squared difference {F n (x) — F(x)} . When ^(x) — 
[F(x){l- F(x)}}- 1 , Eq. © yields the Anderson- 
Darling (AD) statistic. 

The tests derived from the KS and AD statistics are 
called Kolmogorov-Smirnov and Anderson-Darling, re- 
spectively. It is well-known that the KS test exhibits 
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FIG. 3. (Color online) Empirical PDF calculated for simu- 
lated sample from Levy stable distribution with parameters 
a = 1.98, P = 1, a = 1 and n — and for normal distribution 
(with mean equal zero and variance equal 2) in the log-log 
scale. Number of observations in both samples is 5000. The 
straight line shows the asymptotic slope of the Levy stable 
PDF with a = 1.9. Inset: Empirical cumulative fourth mo- 
ment for the samples considered. 

poor sensitivity to deviations from the hypothetical dis- 
tribution that occur in the tails, whereas the AD test 
is considered as one of the most powerful when the fit- 
ted distribution departs from the true distribution in this 
area [28| . Since the fit in the tails is of crucial importance 
in the stable case, AD will be our choice when testing for 
the non-Gaussian stable law. 

For testing the Gaussianity we propose to use the stan- 
dard Jarque-Bera (JB) test [29], which is different from 
the KS and AD. It is a goodness of fit test of whether 
sample data have the skewness and kurtosis matching a 
normal distribution. The JB statistic is defined as: 

where S and K are the sample skewness and kurtosis 
respectively, namely 

■V , "'^ (5) 

k= ^SUO*-*)* Q . (6 ) 
(VV»EJU(^-*) 2 ) 

Equations (H])-©, similarly to Eq. ([1]), can be used for 
any sample from an arbitrary probability distribution. 
The value of the JB statistic given by Eq. (jj) forms a 
random variable which converges to zero if the underly- 
ing distribution has skewness zero and kurtosis 3 (e.g., 



Gaussian). Any deviation from zero skewness and kur- 
tosis equal 3 increases the JB statistic. For distributions 
with infinite kurtosis (e.g., Levy stable with a < 2) it 
diverges to infinity. The test is quite standard and im- 
plemented in various numerical packages, like, e.g., R or 
Matlab. 

To employ any of the tests, first, we need to estimate 
the parameters of the hypothetical distribution. In the 
Gaussian case the standard method is the maximum like- 
lihood, whereas in the non-Gaussian stable case we sug- 
gest the use of the regression method estimates |23l . l3lj . 

The general test of fit is structured as follows. The 
hypothesis is that a specific distribution is acceptable, 
whereas the alternative is that it is not. Small values of 
its test statistic T are evidence in favor of the hypothesis, 
large ones indicate its falsity. To see how unlikely such 
a large outcome would be if the hypothesis is true, we 
calculate the p- value by: p- value = P(T > t), where t 
is the statistic value for a given sample. It is typical to 
reject the hypothesis when a small p-value is obtained, 
like, e.g., below 3% or 5%. To calculate p-values for KS 
and AD tests we apply the procedure proposed in [32| and 
described in detail in (30|. The p- values for the JB test 
were evaluated in Matlab using its standard procedure 
"jbtest". 

Now, we can combine the AD and JB tests with the 
test based on ECFM, to introduce an algorithm for rec- 
ognizing alpha-stable Levy distribution with Levy index 
close to 2. 

In the first step of the algorithm, we propose to ob- 
serve the ECFM of the sample. If the ECFM tends to 
a constant, we check for the Gaussian distribution by 
using the JB test. If its p- value exceeds the confidence 
level (usually 5%), then we can assume the underlying 
distribution of time series is Gaussian. In this case we 
estimate its parameters by using the standard maximum 
likelihood estimation method. If the JB test shows the 
data can not be modeled by a Gaussian distribution, then 
we test them for the Levy stable distribution. If the AD 
test gives the p- value that exceeds the confidence level, 
then we can assume the time series can be descibed by 
the alpha-stable distribution. In this case we estimate its 
parameters via the regression approach (3lj |. 

If the ECFM exhibits chaotic-like behavior, then we 
test for the Levy stable distribution. If the p-value is 
greater than the confidence level, then we can assume 
the data follow the Levy stable law. Its parameters are 
regression method estimates. A scheme of the whole pro- 
cedure is depicted in Fig. |U 

In order to illustrate the procedure we apply it to the 
simulated random samples from Section [XXJ The obtained 
results are presented in Table HI In the last two columns 
we demonstrate the values of the JB and AD statistics, 
respectively, together with the corresponding p- values (in 
parenthesis). 

For the sample from the Levy stable distribution with 
parameters a — 1.98, /3 = 0, a = 1 and p = we can 
observe chaotic behavior of the ECFM, see Fig. [T] There- 
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using the AD test that justifies the stable distribution. 

For the last simulated sample from the normal distri- 
bution we observe that the ECFM tends to a constant. 
The obtained p- value clearly indicates that the underly- 
ing distribution of the analyzed time series is Gaussian. 

The algorithm proposed in the paper provides a pro- 
cedure to distinguish between the case a « 2 and a = 2. 
In order to distinguish the case of a ps 2 and (3—1 with 
a pa 2 and j3 = one has to analyze the estimated skew- 
ness parameter (3. However, we note that the confidence 
interval constructed for this parameter can be quite wide 
if the number of observations is not large enough. This 
circumstance makes the difference not definite. 



IV. APPLICATIONS 



FIG. 4. Schematic algorithm for recognition of Levy stable 
distribution with Levy index close to 2. 



TABLE I. Illustration of the introduced procedure for simu- 
lated samples from Section [TTJ The last two columns contain 
values of the appropriate statistics, calculated by Eqs. ©-([BJl, 
and the corresponding p- values (in parentheses). 



Sample 


ECFM tends 


JB stat. 


AD stat. 


distribution 


to const? 


(p- value) 


(p- value) 


Levy 


NO 




0.24 


a = 1.98, 13 = 






(0.64) 


a = 1, /i = 








Levy 


NO 




0.88 


a = 1.90, 13 = 






(0.42) 


a = 1, fx = 








Levy 


YES 


61.33 


0.18 


a = 1.98, /3 = 1 




(< 0.01) 


(0.86) 


a = 1, fx = 








Gaussian 


YES 


0.55 




fi = 0, a = s/2 




(0.5) 





fore, according to the procedure, we test its stability us- 
ing the AD test statistic. The obtained p- value indicates 
the time series constitutes independent identically dis- 
tributed (i.i.d.) Levy stable random variables. 

As we observe in Fig. [2j the ECFM for the sample 
from symmetric Levy stable distribution with a = 1.9 
and a = 1 behaves chaotically, therefore we use the AD 
statistic here. The obtained p-value clearly indicates the 
underlying stable distribution. 

The ECFM for simulated data from the totally skewed 
stable distribution with a = 1.98 for large number of ob- 
servations tends to a constant (see Fig. [3]), therefore first 
we test the Gaussian law hypothesis. The JB test indi- 
cates the time series can not be modeled by the normal 
distribution, therefore we check for the Levy stable by 



In this section we investigate the data obtained in 
experiment on controlled thermonuclear fusion device. 
Namely, we analyze the statistical properties of plasma 
fluctuations before and after so called L-H transition phe- 
nomenon, that is sudden transition from low confinement 
mode (L mode) to a high confinement mode (H mode) ac- 
companied by suppression of turbulence and rapid drop 
of turbulent transport at the edge of thermonuclear de- 
vice [33j]. The implementation of the H mode regime, 
which is chosen as the operating mode for future ITER 
device [341 ]. requires detailed investigation of the transi- 
tion physics. Electric Langmuir probes rep resent an ef- 
fective diagnostic tool for this purpose [3a|. The impor- 
tant characteristics of edge plasma turbulence, such as 
fluctuation amplitudes, spectra, and turbulence-induced 
transport are investigated in the Uragan-3M (U-3M) stel- 
larator torsatron by the use of high resolution measure- 
ments of density (ion saturation current) and potential 
(floating potential) fluctuations with the help of mov- 
able Langmuir probe arrays [36]. Changes in fluctuation 
behavior are good indicators of transition to improved 
confinement modes. In particular, investigations of spon- 
taneous LH-likc transition (hereinafter, transition) in the 
U-3M torsatron reveal a distinct decrease of fluctuation 
level and radial turbulent particle flux [37| . In this sec- 
tion we analyze the three data sets. They present typi- 
cal floating potential fluctuations (in volts) in turbulent 
plasma, registered in the U-3M torsatron for the small 
torus radial positions r — 9.9 cm (data set SI), r = 10 
cm (S2), and r = 11.25 cm (S3). For the detailed descrip- 
tion of the experimental set-up see Refs. [H], HH- One 
can observe obvious change of the fluctuation amplitude 
in the middle position of every data source, which indi- 
cates the transition [37]. Therefore, we split each data 
set into two parts: before (datal) and after the tran- 
sition point (data2). Moreover, we examine only such 
parts of each data which constitute the stationary pro- 
cess. Precisely, each datal contains observations between 
6000 and 12000, whereas data2 between 14000 and 19000. 
ADC sampling rate is 625 kHz that means the time res- 
olution between two consecutive measurement points is 
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1.6 mcsec. In Fig. [5] we depict the examined data sets 
and the selected parts before and after the L-H point. 
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FIG. 5. (Color online) The empirical time series from plasma 
physics described in Section HVl 

According to the procedure of testing between Levy 
stable and Gaussian distributions presented in Section 
IIII1 we start with examining the ECFMs for all data sets. 
The result of this simple test together with right tails of 



the empirical PDFs are presented in Fig. [6] 
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FIG. 6. (Color online) The right tail of the empirical PDFs 
and ECFMs for six examined data sets. 

This visual test indicates that following data sets can 
not be considered as a Gaussian time series: SI - datal, 
S2 - datal and S3 - datal, whereas for the rest of the 
examined samples, i.e., SI - data2, S2 - data2 and S3 
- data2, the ECFM stabilizes. In order to confirm this 
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TABLE II. Value of the test statistic and corresponding p- 
values (in parentheses) for the Gaussian law hypothesis (the 
JB test). 



Data set 


JB stat. 




(p- value) 


SI - datal 


12.35 




(< 0.01) 


SI - data2 


7.39 




(0.025) 


S2 - datal 


71.59 




(< 0.01) 


S2 - data2 


1.67 




(0.43) 


S3 - datal 


62.25 




(< 0.01) 


S3 - data2 


0.96 




(0.5) 



fact, we test all examined data sets for the Gaussian dis- 
tribution by using the JB test presented in Section IIIII 
The p- values and values of the JB statistics are depicted 
in Table HU 

The JB test clearly indicates that SI - datal, S2 - 
datal and S3 - datal follow a non-Gaussian distribution 
because the p- values are smaller that the confidence level 
5% (in fact they are even less than 1%). Let us now con- 
sider the sample SI - data2. In this case the ECFM tends 
to a constant whereas the JB test rejects the hypothesis 
of Gaussianity with p-value close to the confidence level. 
This is why we use another test in order to confirm one 
of the hyphotheses (Gaussian or non- Gaussian distribu- 
tion). We propose to use the AD test presented in Sec- 
tion IIIII for testing for Gaussianity. The value of the 
AD statistic is 0.5125, while the corresponding p- value is 
equal to 0.16, which indicates the Gaussian behavior. 

In the next step of our analysis we test the samples 
with chaotic ECFMs for the Levy distribution, namely 
SI - datal, S2 - datal and S3 - datal. In Table HIT1 we 
present p- values and values of the AD statistics for the 
hypothesis of the stable distribution that indicate stable 
behavior of examined data sets. 

Finally, we can present the index of stability and the 
skewness parameter for samples SI - datal, S2 - datal 
and S3 - datal. They were obtained by the regression 
method (23|. The estimated (3 parameter for SI - datal 
is equal to 0.2 while for S2 - datal is on the level 1 
and for S3 - datal is —1. The results of the estimated 
parameter a and the corresponding confidence intervals 
are presented in Table IIV1 



TABLE III. Values of the test statistics and corresponding 
p-values (in parentheses) for the Levy stable law hypothesis 
(AD test). 



Data set 


AD stat. 


SI - datal 


0.39 
(0.24) 


S2 - datal 


0.48 
(0.72) 


S3 - datal 


0.62 
(0.91) 



TABLE IV. Estimated parameter a of the Levy stable distri- 
bution with the corresponding confidence intervals. 



Data set 


a 


SI - datal 


1.98 
[1.97; 1.99] 


S2 - datal 


1.97 
[1.96; 1.99] 


S3 - datal 


1.98 
[1.96; 1.99] 



V. DISCUSSION 

Let us first add a word of caution: of course, as with 
any other family of distributions it is not possible to prove 
that a given data set is or is not stable [38j . and even 
testing for normality is still an active field of research [3!| . 
In the empirical data analysis like the one presented here 
the best what can be done is to determine whether or not 
the data are consistent with the hypothesis of stability. 

However, with the testing procedure developed in 
present paper we were able to distinguish between Gaus- 
sian and stable distribution with Levy index close to 2. In 
other words, we showed that the difference between such 
laws is essential and one can identify the alpha-stable dis- 
tribution with a close to 2, even if the number of available 
data points is not too large. We illustrated the procedure 
on simulated samples from Gaussian and non-Gaussian 
stable distributions. The results show that the algorithm 
is reliable and works well. It might be useful for detecting 
Levy stable distributions in different experiments. 

As an example, we applied the developed procedure 
to assess the stable properties of the fluctuations mea- 
sured in the edge plasma of stellarator torsatron "Uragan 
3M" and demonstrated that the statistics of fluctuations 
changes from stable to normal at the LH transition. 

Another remark concerns the observed phenomenon of 
changing statistics. We first note that the Levy statistics 
of plasma fluctuations has been detected in the measure- 
ments on stellarators "Uragan 3M" and Heliotron J by 
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using the method of quantiles [T3| and in the tokamak 
ADITYA by measuring probability to stay at the origin 
[lij . The Levy indexes of the detected stable distribu- 
tions were sufficiently less than 2. Non-Gaussian heavy 
tailed PDFs for low frequency turbulence have been also 
observed in other toroidal plasma confinement systems 
such as T-10 tokamak, L-2M, TJ II and LHD stellarators 
(ioj as well as in astrophysical plasmas [U [42[. As to 
the physics of the change of the statistics, it in turn, im- 
plies qualitative changes of the basic properties of plasma 
turbulence, and should be taken into account in the theo- 
retical models of the LH transition. This important issue 



is out of the scope of the present paper. On the other 
hand, we believe that the change of statistics that we ob- 
served in the data from "Uragan 3M" may inspire other 
authors to check whether LH transition is accompanied 
by change of statistics in other plasma confinement de- 
vices as well. 
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